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Neutron stars are efficient accelerators for bringing charges up to relativistic energies. We show 
that if positive ions are accelerated to ~ 1 PeV near the surface of a young neutron star (t age ;$ 10 s 
yr), protons interacting with the star's radiation field will produce beamed /i neutrinos with energies 
of ~ 50 TeV that could produce the brightest neutrino sources at these energies yet proposed. These 
neutrinos would be coincident with the radio beam, so that if the star is detected as a radio pulsar, 
the neutrino beam will sweep the Earth; the star would be a "neutrino pulsar". Looking for 
emission from young neutron stars will provide a valuable probe of the energetics of the neutron 
star magnetosphere. 

PACS numbers: 97.60.Gb, 94.30.-d, 95.85.Ry 
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Pulsar emission is interpreted as a consequence of the 
acceleration of charges to relativistic energies somewhere 
in the neutron star's magnetosphere. As the charges 
(positive ions or electrons) move along the curved mag- 
netic field lines, they produce curvature radiation which 
in turn produces a pair cascade and beamed radiation 
(see for a review). A long-standing and unsolved 
problem in neutron star magnetospheric physics is where 
charges obtain most of their energy. Near the stellar 
surface, charges are constrained to move along magnetic 
field lines. In the presence of the plasma, the field lines 
are approximate equipotentials; the electric and magnetic 
fields satisfy E • B ~ for a quasistatic magnetosphere. 
In most theories of pulsar emission, charge acceleration 
occurs in a charge depleted polar region very near the 
stellar surface, sustained by a deviation of the magneto- 
sphere from corotation with the star 0- C3- E3| ■ An alter- 
native is the "outer gap" model [j| in which the acceler- 
ation site can be far from the stellar surface. To deter- 
mine exactly where and how the charges are accelerated 
would require a detailed description of the neutron star 
magnetosphere, which we currently lack. Here we show 
that young neutron stars might be among the brightest 
sources of muon neutrinos yet proposed, and so might be 
the first sources to be detected by planned or operating 
neutrino telescopes such as AMANDA-II and IceCube 
0, and ANTARES, NEMO and NESTOR 0. The de- 
tection of young neutron stars as strong neutrino sources 
would be a fascinating breakthrough in its own right, and 
would indicate that protons (or heavier ions) are acceler- 
ated to high energies close to the star, thus providing an 
important probe of the energetics and physical conditions 
that prevail in the poorly understood magnetosphere of 
a neutron star. 

If the star's magnetic moment fj, and angular velocity 
n satisfy fi ■ £1 < (as one would expect for half of 
neutron stars), positive ions will be accelerated to infinity 



If the neutron star is young, its surface will emit 
brightly in soft x-rays, and the protons in accelerated 
nuclei will scatter with this radiation field. If the protons 
are sufficiently energetic, they will exceed the threshold 
for photomeson production through the A resonance (the 
A + is an excited state of the proton, with a mass of 
1232 MeV). The A+ quickly decays to a tt+, and muon 
neutrinos are produced through the following channels: 
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This idea for production was explored by |9( in the 
context of magnetars. In this paper we explore a similar 
idea in a different physical regime, that of young, hot 
neutron stars. We show that young neutron stars could 
be strong sources of muon neutrinos, with energies ~ 50 
TeV, and with fluxes observable by large-area neutrino 
observatories. 

Energetics. Goldreich & Julian [j| showed that the 
potential drop across the field lines of a pulsar from the 
magnetic pole to the last field line that opens to infinity 
is of magnitude A$ = /iil 2 /2c 2 , where fx — BR 3 is the 
stellar magnetic moment, B is the strength of the dipole 
component of the field at the magnetic poles and R is the 
stellar radius. If charges in rapidly rotating pulsars with 
unexceptional fields (B ~ 10 12 G) could be accelerated 
by this potential drop, the corresponding energy would 
be huge (neglecting radiation losses, but see below) 



e$ = eZA$ = 7ZB 12 p~l EeV, 



(2) 



per ion, where e is the electronic charge, Z is the charge 
number of the ion, Bn = B/10 12 , p ms is the spin pe- 
riod in milliseconds and R ~ 10 6 cm is the stellar radius. 
Near the star, charges are constrained to follow the mag- 
netic field lines. For a charge to be accelerated, E • B / 
is required. Moreover, there must be at least some ac- 
celeration near the stellar surface, to supply the charges 
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that support the magnetosphere. As a conjecture to be 
tested by neutrino observations, we shall assume that a 
strong, accelerating field exists near the stellar surface. 
In this case, radiation losses will determine whether or 
not nuclei can attain energies of the magnitude given by 
cq. [2| . The power dissipated by a charge q of energy 
e moving along a field line of radius of curvature p c is 
de/dt — (2/3)(q 2 c/ p 2 )(e/mc 2 ) A . The acceleration time 
scale is ~ p c /c, giving de/dt ~ ec/p c . Near a pole of a 
magnetic dipole, p c > R. Radiation losses limit the ac- 
celeration per nucleon of a nucleus of mass number A 
to be e max ~ 20(A/Z 2 ) 1 / 3 PeV. Hereafter we assume 
Z/A~ 1/2. 

Now we estimate the photomeson production threshold 
to compare with the radiation loss limit. The threshold 
condition for a proton to reach the A resonance is that 
the proton and photon energies, e p and e 7 , respectively, 
satisfy 

e p e 7 > 0.3 GeV 2 / 3 , fg = 0-~ cos^)" 1 , (3) 

where 9 pl is the incidence angle between the proton and 
the photon in the lab frame. This condition is indepen- 
dent of whether the proton is free or bound in a nucleus; 
in either case, photomeson production occurs in nearly 
the same way. In the rest frame of the nucleus, the pho- 
tons have an energy of about 300 MeV. This energy is 
much larger than the binding energy per nucleon, so nu- 
clear binding has little effect on the conversion of protons. 

Young neutron stars (t age < 10 5 yr) typically have tem- 
peratures of ^ 0.1 keV (= 1.2 x 10 6 K). Typical pho- 
ton energies near the surface are e 7 = 2.8fcToc(l + z g ) ~ 
0.4 keV, where z g ~ 0.4 is the gravitational redshift. The 
proton threshold energy for the delta resonance is then 
e„,th - T-L v f g PeV, where T .i kcV = (fcT^/0.1 keV). 
As we show below, the conversion, if it happens, must oc- 
cur near the stellar surface. Hence, protons scatter with 
photons with a lab frame scattering angle of 9 P1 ~ 90° , 
giving f g ~ 1. (We neglect the effects of gravitational 
light bending.) In comparing the threshold energy e p ,th 
to the radiation loss limit e max , it appears that radia- 
tion losses will not preclude protons from reaching the A 
resonance, even for heavy nuclei. 

Is the accelerating potential strong enough to bring 
protons to the A resonance? The maximum energy per 
proton is e$/A. Since the protons are scattering with 
radiation that is being radiated isotropically from the 
star (we ignore the possibility of temperature variations 
across the surface), 9 P1 is < 90°, and so f g > 1. For the 
A resonance to be reached, we require that 

Bi 2 p^T 0AkeV > 3 x 10~ 4 . (4) 

This condition, provided that /j, ■ ft < 0, represents the 
most optimistic condition for neutrino production since 
we have assumed that the full potential A<i> is available 
for accelerating ions. Eq. 0] is satisfied in many young 



pulsars, though not all. Magnetars, which have typical 
fields of 10 15 G, spin rates of 4 — 8 s and T .ik o v — 4 — 6, 
are borderline neutrino emitters according to eq. [2[, as 
found by 0. In these estimates, we have neglected the 
possible quenching of the accelerating field by pair cas- 
cades produced through inverse-Compton scattering of 
the protons with the radiation field; numerical simula- 
tions of this effect indicate that the quenching will be 
negligible, even for the high surface temperatures we are 
considering |l0j |. Failure to detect neutrinos from pulsars, 
however, could mean that the field is strongly quenched. 

Conversion Probability. For a significant neutrino 
flux to exist, the conversion probability for pj — » A + 
must be sufficiently high. The value of f g averaged over 
the surface increases rapidly with distance from the star, 
hence, the conversion must occur near the stellar surface 
or the threshold becomes unattainable. We suppose that 
the A resonance is reached before an altitude at which 
photons with f g <2, become unavailable. The photons 
with the lowest values of f g originate from the stellar 
limb. To crudely estimate f g , we consider 9 P1 to be de- 
termined by the angle between the center of the star and 
the stellar horizon as "seen" by a proton at r > R, mea- 
sured from the stellar center. With this definition, f g = 2 
corresponds to r = 1.2i?. Some models of charge accel- 
eration in pulsar magnetospheres find that significantly 
charge-starved regions can exist just above the stellar 
surface (~ 10 4 cm), so that the charges receive most of 
their acceleration there 0- 0- E| ■ well below the altitude 
at which f g begins to exceed 2. 

The mean free path for conversion at radius r is 
l(r) = (n 7 cr P7 ) • , where n 7 is the number density of 
the radiation field at radius r and a pl ~ 5 x 10~ 28 cm 2 is 
the cross section for A + production. Assuming isotropic 
radiation from the stellar surface, the photon density 
there is u 7 (i?) = (a/2.8fc)([l + z^T^) 3 ~ 9 x 10 19 T 3 lkcV 
cm~ 3 , where a is the radiation density constant, and 
n 7 (r) = n 1 {R)(R/r) 2 . The mean free path at the stellar 
surface is l(R) ~ 2 x 10 7 T~ 3 kcV cm. The probability P 
that a proton has not been converted changes with r as 
dP/P = —dr/l(r). By the time a proton starting at the 
surface has reached radius r, the probability of conver- 
sion is Pconv(^) = 1 — P{ r )- Requiring conversion in the 
range R <r < 1.2R gives P conv ~ 0.02T 3 lkcV . For the 
surface temperatures we are interested in (T .ik o v — 1), 
the radiation field is quite optically thin to the conver- 
sion process. Hence we can expect the conversion of only 
one proton per nucleus, if there is a conversion at all. Be- 
cause the energy imparted to the A + is much larger than 
the binding energy per nucleon, the A + will be ejected 
from the nucleus. The pulsar emission mechanism will 
be essentially unaffected, since so few nuclei are affected. 
The probability of photodisintegration is negligible. 

The Neutrino Energy. The average fraction of en- 
ergy transferred from the proton to the pion is <~ 0.2, 
or 200 Tq"^ v TeV, corresponding to a Lorentz factor 
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of 7?r = 1.4 x 10 6 Tjy lkeV . The pion decays in a time 
7"d = t^tt, where 7V = 2.6 x 10~ 8 s, and thus moves to 
r ~ 10 3 i? before decaying. The pion is in the densest 
part of the radiation field for only a short time, of or- 
der R/c ~ 3 x 10~ 5 s, and suffers negligible energy loss 
to inverse-Compton scattering. When the pion decays, 
it gives one-fourth of its energy to the muon neutrino, 
the rest going to the other three leptons. The resulting 
neutrino energy is then 

e„„ * 50 T-^ v TeV. (5) 

The Neutrino Flux. The accelerated protons are 
far more energetic than the radiation field with which 
they are interacting. Any pions produced through the 
A resonance, and hence, any muon neutrinos, will be 
moving in nearly the same direction as the proton was 
when it was converted. We thus assume that the star 
produces a neutrino beam that is nearly coincident with 
the radio beam, that is, if the star is detected as a radio 
pulsar, the neutrino beam will sweep the Earth. The 
radio pulsar will also be a "neutrino pulsar". 

For a quasistatic magnetosphere, the charge density 
near the stellar surface is p q ~ eZn — B/pc (cgs), where 
no is the Goldreich- Julian number density of ions ||. 
To get acceleration of charges, there must be a charge 
depleted gap somewhere above the star. We parameterize 
the density in the gap as fdJiQi where fd < 1 is a model- 
dependent and uncertain depletion factor. The depletion 
factor should be regarded as a free parameter, though 
even modest depletion (f d ~ 1/2) could give enormous 
accelerating fields, and a significant fraction of e$ could 
be reached 0, 0- We see the radio beam for only a 
fraction fb of the pulse period, the duty cycle. Typically, 
ft, ~ 0.1 — 0.3 for younger pulsars. We take the duty 
cycle of the neutrino beam to be fb- The phase-averaged 
neutrino flux at Earth resulting from the acceleration of 
positive ions, at a distance d from the source, will be 

$v - Cfbfd.no (J^j -Pconv (6) 

Large-area neutrino detectors use the Earth as a medium 
for conversion of a muon neutrino to a muon, which 
then produces Cerenkov light in the detector. The 
conversion probability in the Earth is P v ~ 1.3 x 
10 _6 (e„ ll TeV) where e„ is the energy of the incident 
neutrino Combining eqs. [S] and [Hj gives a muon 

event rate of 

dN 

— - 4>uPu^ 

~ 10 5 Z- l f b f d B 12 p-ld^ c Tl lkeY km- 2 yr- 1 . (7) 

As an example, if T .i k cV = !,/& = 0.1, fd = 0.1, p ms 
1 0. d kpc = 3, Z = 2, the muon flux would be ~ 6 km -2 
yr -1 at a neutrino energy of ~ 50 TeV. The rate given by 



Source 




age 


Pms 


B\2 


Tb.ikcv 


f b dN/dAdt 






yr 








km - yr -1 


Crab 


2 


10 3 


33 


3.8 


< 1.7 [13] 0.14 1200 


Vela 


0.29 


io 4 - 2 


89 


3.4 


0.6 [14] 


0.05 800 


J0205+64 3.2 


1Q 2.9 


65 


3.8 


< 0.9 [15] 


0.04 20 


B1509-58 


4.4 


1Q 3.2 


151 


15 


1? 


0.26 130 


B 1706-44 


1.8 


IO 43 


102 


3.1 


1? 


0.13 120 


B1823-13 


4.1 


1Q 4.3 


101 


2.8 


1? 


0.34 60 


Cass A 


3.5 [16] 


300 


10? 


1? 


4 [16] 


0.1? 1300 


SN 1987a 


50 


17 


1? 


1? 


4? 


0.1? 60 



TABLE I: Estimated fi Fluxes at Earth. Magnetic fields 
were estimated from the vacuum dipole prediction: B = 
(3.2 x 10 39 pp) 1//2 G. Numbers followed by question marks in- 
dicate guesses. We took fd — 1 and Z = 1 so these rates are 
upper limits, and must be scaled by a factor fdZ -1 . Radio 
pulsar spin parameters were taken from [r?ll . Temperatures 
and limits on temperatures were taken from the references 
indicated. The temperature upper limits on the Crab and 
J0205+64 were used. 

eq. is model independent and quite general, provided 
that young neutron stars succeed in accelerating positive 
ions to the A resonance. 

Detectability. For a neutron star to be a detectable 
neutrino pulsar, it must satisfy eq. [4], and be sufficiently 
close to give a strong flux. For a neutron star to be a de- 
tectable neutrino pulsar, it must satisfy -Bi2Pms^o.ikcV > 
3 x 10~ 4 (eq. 0]), and be sufficiently close as to give a 
strong flux. Potentially interesting sources and upper 
limits for their muon event rates are given in Table 1. 
We have optimistically included SN1987a in the Table, 
in case this supernova remnant contains a neutron star. 
We have estimated the background above 50 TeV due 
to atmospheric cosmic ray events using the flux given by 
[13 ; for detectors such as AMANDA-11 and ANTARES, 
with angular resolutions of ~ 1°, the atmospheric back- 
ground above 50 TeV contributed over several years will 
be so low that the detection of a single muon would rep- 
resent a statistically significant neutrino detection. 

Discussion. Results of 607 d of data from AMANDA- 
II are now available [lilllfij . Cass A and the Crab pulsar 
were not detected as neutrino sources. Cass A has not 
been observed as a pulsar, and so would not be expected 
to be a source in this model. The Crab pulsar would be 
a strong source, but only if fi ■ ft < 0. If the Crab is a 
neutrino pulsar, and the temperature is close to the ob- 
servational upper limit, lack of detection by AMANDA-II 
requires /d;$0.01 [Z = 2), which seems improbably low 
in the face of current models of charge depleted gaps. 
It seems likely that the Crab is simply not a neutrino 
pulsar, perhaps because n ■ fl > 0. (The threshold con- 
dition eq. 0] would be easily met even for much lower 
temperatures than the upper limit given in Table 1). 

Should AMANDA-II have seen neutrinos from other 
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pulsars? There are now 64 known pulsars that are within 
10 kpc of Earth and younger than 10 5 yr 0, detected 
mostly by the Parkes Radiopulsar Survey. Any of these 
is a potential neutrino pulsar. Because the Parkes tele- 
scopes are located in the southern hemisphere, only 13 
of the 64 radio pulsars are in the northern sky, the 
part visible to AMANDA-II. To make an estimate of the 
probability of detection of a random neutrino pulsar by 
AMANDA-II, suppose that AMANDA-II were located at 
the North Pole, observing southern-sky sources. We can 
then use a sample of the 51 pulsars from the southern sky 
to assemble a probability distribution for the rates. We 
took To.ikcV = 1, fd — 0.5 and Z = 2 in all cases, /& in 
the cases that it has been measured, and /& = 0.1 in all 
other cases. Using these uncertain, but very conservative 
numbers, we find only ~ 1 pulsar that satisfies eq. 0] and 

■ O < 0, so it is not surprising that AMANDA-II has 
not yet seen a neutrino pulsar. Given the small number 
of pulsars known in the northern sky, our model cannot 
be statistically constrained by AMANDA-II. AMANDA- 
II's successor, IceCube, with ten times the effective area, 
might reveal neutrino pulsars. Moreover, in the south- 
ern sky there are a number of promising candidates 
that might be detected by ANTARES (expected to go 
into operation during 2005), but that are inaccessible to 
AMANDA-II; these include Vela, B1509-58, B1706-44, 
B1823-13 and possibly SN 1987a. There are other candi- 
dates in the southern sky, none with measured temper- 
atures, but assuming To.ikcV > U there are nine known 
pulsars that satisfy eq. @], are within a distance of 5 
kpc and are younger than 10 5 yr and therefore likely to 
have To.ikcV ~ 1- These give typical muon event rates of 
~ fdZ~ 1 (20 — 130) km -2 yr -1 ; three of these candidates 
are listed in Table 1. 

The most promising candidate source in the southern 
sky is the Vela pulsar, which if detected by ANTARES, 
would give an event rate of 80Z~ 1 fd yr -1 . For Z = 2, 
the source would be detectable over a year of observing 
time even if the acceleration region is charge depleted to 
~ 3% of the Goldreich- Julian density. For modest deple- 
tion (fd ~ 1/2), Vela could be detected in under three 
weeks. If we are unfortunate, and Vela has /Lt • CI > 0, it 
will not be a neutrino pulsar, and we must look for weaker 
sources. For the other candidates to be detectable by 
ANTARES, B 1508-58 and B 1706-44 might be detected 
in about four months for fd — 1/2 and Z = 2. After 
the first year of data from ANTARES, statistical con- 
straints will be possible even if there arc no detections, 
since the the properties of pulsars in the southern sky 
are well understood. No detections would mean that 
photomeson production is negligible or nonexistent, pro- 
viding an interesting constraint on the energetics of the 
magnctospherc. A detection would allow a determination 
of Z~ 1 fd, and hence, the charge density in the accelera- 
tion region. We stress that the rates appearing in Table 
1 are upper limits; they could be significantly reduced by 



the factor Z^ 1 fd- 

We conclude by pointing out that neutrino pulsars 
could be among the strongest sources of high-energy 
muon neutrinos in the sky, and therefore, might be the 
first sources to be detected. (See [2(3, for a review of 
other models). We will present the energy spectrum of 
neutrino pulsars in a forthcoming publication. 
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